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NOTATION 


a semimajor axis, km 

e eccentricity 

H altitude, km 

L minimum angle between departure vector and orbit plane, deg 
i inclination, deg 

J 2 second harmonic of planetary oblateness 

R planet radius, km 

t time, sec 

AV propulsive velocity increment, km/sec 

V„ hyperbolic excess velocity, km/sec 

x,y,z components of the planet centered coordinate system 
6 declination, deg 

r unperturbed orbital period, sec 

p right ascension, deg 

£2 longitude of ascending node, deg 

co argument of periapsis, deg 

Subscripts 

h hyperbola 

i intermediate or departure ellipse 

o original parking orbit ellipse 

p periapsis 
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ANALYSIS OF LAUNCH WINDOWS FROM ELLIPTICAL ORBITS 


FOR REPRESENTATIVE MARS MISSIONS 1 

Larry A. Manning and Byron L. Swenson 

Office of Advanced Research and Technology 
Mission Analysis Division 
Moffett Field, California 94035 


SUMMARY 


Round trips to Mars have been investigated to define representative launch windows (i.e., 
allowable staytime variations) for Mars departures from elliptical orbits and the associated 
propulsive velocity requirements. The 1982 inbound and 1986 outbound Venus swingby 
opportunities were selected for analysis and serve to demonstrate the influence of the 
interplanetary trajectory characteristics on the launch window propulsive velocity requirements. 
These velocity requirements are the sum of the Mars arrival and departure maneuvers. This report 
presents contour maps indicating the effects on the launch window of velocity capability 
(total AV), transfer technique, orbit inclination, and orbit eccentricity. Use of one-, two-, and 
three-impulse transfers from the Mars orbit to the Earth return trajectory were investigated. For all 
cases, insertion at planet arrival was into an orbit coplanar with the arrival asymptote and any 
required plane change was performed during the planet departure phase. 

Two specific conclusions can be drawn from the results of this study. First, as a result of the 
effectiveness of the three-impulse transfer in reducing the AV requirement to a nearly constant 
value regardless of the orbit geometry, the orbit selection can be optimized from other 
considerations (i.e., science, communication, power, etc.) independent of the interplanetary 
trajectory considerations. Second, for highly elliptical orbits (e = 0.7) a total AV (arrival plus 
departure) at Mars that is 1 km/sec below the minimum coplanar requirement for circular parking 
orbits provides staytimes of 60 days and usually allows any orbit inclination to be used. This AV is 
effectively 2 km/sec below that required for circular orbits with the same allowable staytime. In 
addition, whereas the circular orbits require three impulses at departure, for some missions a single 
impulse can be used to eject from the elliptical orbit. 


INTRODUCTION 


Elliptical parking orbits are receiving increased attention for use as planetary parking orbits 
during interplanetary exploration. These orbits have been shown to greatly reduce the propulsive 
requirements for arrival and departure at the target planet (ref. 1 ). In addition, orbital plane 

1 The material in this report was summarized in a paper by the same authors and Jerry M. Deerwester, and 
entitled “Launch Window Analysis for Round Trip Mars Missions,” presented at the Canaveral Council of Technical 
Societies’ Fifth Space Congress, Cocoa Beach, Florida, March 11-14, 1968. 



changes can be performed with a lower propulsion system requirement for elliptical orbits than for 
the low circular orbits generally considered for planetary exploration. 

Selection of the orbit inclination at the planet is made after consideration of the mission 
science requirements (optical coverage, probes, landers, etc.) and of the effect of that inclination on 
the mission propulsive requirements. This second effect is particularly significant since the orbit 
plane will not, in general, contain both the arrival and the departure velocity vectors. Thus it is 
necessary to perform a plane change maneuver, which can be costly in terms of propulsive 
requirements. 

This report presents the results of an analysis of the influence of elliptic orbit inclinations on 
the staytime in orbit and the associated energy requirements for representative round trip Mars 
missions in a similar manner previously used for circular orbits (ref. 2). In addition to the orbit 
itself, the technique used to transfer between the orbit and the interplanetary trajectory is 
analyzed. Three techniques are considered. They are: (1) single impulse; (2) two impulse, the first 
of which provides all the plane change and the second all the velocity change, and (3) three impulse, 
which uses an intermediate ellipse of high eccentricity from which all the plane change is made. 

To ensure realistic results, specific missions were selected for the analysis. The strong coupling 
of the orbit and both legs of the interplanetary transfer is such that an analysis of the orbit transfer 
technique by itself, while interesting, does not provide the insight necessary to allow proper 
mission/spacecraft tradeoffs. Two missions, encompassing the four trajectory types of primary 
interest (i.e., direct outbound and inbound, and Venus swingby outbound and inbound), were 
chosen as representative of future round trips of interest. While the numerical results obtained are 
not directly applicable to other Mars missions, the data trends and resulting conclusions are 
generally applicable. 

The report is divided into two major sections. The first describes the basic parts of the study, 
that is, the orbit geometry, the orbit transfer techniques, and the nominal missions descriptions. 
The second couples the representative missions with the orbit characteristics to provide contour 
maps of constant velocity increments as functions of inclination and staytime in orbit. 


METHOD OF ANALYSIS 


Orbit Geometry 

In considering the launch window problem for a given mission, it is first necessary to establish 
the relative positions of the orbit plane and the hyperbolic escape asymptote with time. For the 
analysis of planetary launch windows, it was assumed that the parking orbit at planet arrival was 
coplanar with the arrival asymptote and was established by a tangential impulse not necessarily at 
periapsis. The resulting orbit elements are shown in figure 1. As indicated, a planet-centered 
right-hand coordinate system with Z axis at the north pole and X axis at the planet vernal equinox 
was chosen. The arrival asymptote or hyperbolic excess velocity vector is defined conventionally as 
a planet centered vector with right ascension p and declination 5. The orbit elements of interest, 
that is, the inclination i and longitude of the ascending node £2 are related by 


2 



sin(p - Sf) 


tan 6 
tan i 


for all i > 5, while the argument of periapsis co is measured in the orbit plane from the ascending 
node. 



Figure 1.- Orbit geometry at planet arrival. 

It can be seen that for each inclination there are two orbit planes coplanar with the arrival 
asymptote. These two orbits can be distinguished by the relative positions of the spacecraft 
approach vector at the time of the final midcourse maneuver and the planet centered excess 
velocity vector. If the spacecraft approaches above the excess velocity vector, it moves initially 
toward the north pole of the planet and the magnitude of (p - £2) will be greater than 90°. This 
orientation will be referred to as a northern insertion. In the other case, the spacecraft approaches 
below the excess velocity vector (i.e., toward the south pole) and the magnitude of (p - £2) is less 
than 90°. Figure 1 illustrates this configuration, which is called a southern insertion. 



Figure 2.- Orbit geometry at planet departure. 


The situation at some time after arrival 
is illustrated in figure 2. The perturbation due 
to planet oblateness causes the orbit plane to 
regress about the planet in the manner shown. 
That is, the inclination of the orbit remains 
unchanged, the longitude of the ascending 
node changes by A£2, and the argument of 
periapsis changes by Aoj where, to first 
order in the planet oblateness (ref. 3) 

A£3 = 37rJ 2 ^j^ (cos i)N 

2 

Aw = (2 - sin 2 N 


3 


with the number of perturbed orbits (N) and the semilatus rectum (p) defined by 

2 


N = 


i 3 T 
1 ♦ 2 J 2 


(I) 0 - 1 sin2 o » 


e 2 ) I/2 


a(l - e 2 ) 


At that time, the orbit plane makes an angle with the departure asymptote. This is the angle 
that must be compensated for during the departure manuever. 


Transfer Techniques 

Three transfer techniques were considered for the departure maneuver: an optimum single 
impulse, a constrained two impulse, and a three-impulse transfer. 

One-impulse technique— In this technique, since at a given time both the departure hyperbolic 
excess vector and the characteristics of the parking orbit at departure are fixed, the optimum 
noncoplanar, nontangential departure maneuver true anomaly is established by a numerical method. 
In general, the minimum total AV for arrival plus departure requires an off-periapsis arrival and 
both an off-periapsis and a nonminimum plane change angle maneuver at departure. The associated 
equations are described in appendix A. 

Two-impulse technique- The first impulse is applied at apoapsis to perform the required 
plane change. The second impulse is applied tangentially near periapsis and provides the required 
velocity change to achieve the desired excess velocity. The associated equations are described in 
appendix B. 

Three-impulse technique- The first impulse is applied at periapsis to increase the orbit 
eccentricity; a value of 0.9 was used for the eccentricity of the intermediate ellipse. The second 
impulse performs the necessary plane-change maneuver at apoapsis, and the third impulse applied 
near periapsis provides the additional velocity increment necessary to achieve the desired hyperbolic 
trajectory. The second and third impulses of this technique are computationally the same as those 
in appendix B for the two-impulse transfer. The first impulse serves to raise the apoapsis altitude, 
thus reducing the subsequent plane change velocity increment for a given 1^. This sequence is an 
approximation of the optimal velocity/plane change distribution between the three maneuvers 
(ref. 4). 


Nominal Missions 

To obtain meaningful operation data and conclusions on launch windows, it is necessary to 
define reasonably representative round trip missions. The missions selected as representative of 
those most likely to be used are the 1982 inbound Venus swingby and the 1986 outbound Venus 
swingby. The study assumed a single outbound (Earth to Mars) leg for each mission independent of 
the staytime. The inbound (Mars to Earth) leg was defined as a function of staytime for each 
opportunity considered, to establish the Mars departure velocity vector variation. The mission 
selection procedures and criteria are discussed fully in a previous report on circular orbits (ref. 2). 
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The outbound leg characteristics are delineated in table 1 , while the characteristics of the inbound 
legs are shown as a function of staytime in figures 3(a) and (b). The break in the curves of 
figure 3(b) occurs where the central angle of the 1986 inbound leg is nearly equal to 180°. 


TABLE 1.- NOMINAL OUTBOUND LEG CHARACTERISTICS 


Leave Earth 

V*, I AV , b 



p, 

deg 

<5, 

deg 

145.4 

74.5 

-0.6 

19.9 


Pass Venus 


P 9 


Arrive Mars 


Date 3 - 

Vco 

P 

5210 

6500 

3.7 

6.0 

288.1 

285.0 



All dates are Julian minus 244 0000 
150-km circular orbit 


Right ascension 


^ 280 — o -20 — £260 




40 60 

Staytime, days 


(a) 1982 inbound swingby (Mars arrival date = 244 5210). 
Figure 3.- Inbound leg trajectory characteristics. 













(b) 1986 outbound swingby (Mars arrival date = 244 6500). 
Figure 3.- Concluded. 


Impulsive Maneuvers 

This analysis assumes that the velocity increments are applied impulsively, and one effect of 
this assumption is the neglect of gravity losses. Gravity losses, however, are small (i.e., 5 percent of 
the total velocity requirement for an I S p of 1000 sec (ref. 5) decreasing to about 3 percent for a 
chemical system I S p of 400 sec (ref. 6)) and can be neglected without significant error. It has also 
been shown that the orientation between the parking orbit and the escape hyperbola is preserved 
whether the velocity addition is impulsive or through a finite thrusting time, although the true 
anomaly of ignition point varies slightly between the two assumptions (ref. 2). 


6 



NUMERICAL ANALYSIS 


Launch Window Development 

In this study, the sum of the arrival and the departure AV requirements was used as a 
measure of the relative performance of the three transfer techniques and to define the staytime 
variation (launch window) for a given energy level. This summation approach is necessary since the 
location of the elliptical orbit periapsis relative to the arrival or departure velocity vector exercises a 
major influence on the magnitude of the transfer AV. 

A single-impulse coplanar insertion maneuver was assumed at arrival. If this maneuver is a 
tangential, periapsis insertion (i.e., minimum arrival AV), the departure maneuver, after a given 
staytime, may require an excessively large AV to compensate for its being performed significantly 
off the optimum departure location. An off-pericenter insertion maneuver at arrival can reduce the 
off-pericenter requirement at departure and thereby lower the sum of the AV. This is a reflection 
of the fact that coplanar maneuvers as much as 30° to 40° off-periapsis can be made for 
a AV penalty of about 2 percent or less (ref. 2). The value of the optimal off-pericenter angle at 
arrival is a function of the inclination of the parking orbit and the desired staytime at the planet. 

To assess the launch window requirements, a nominal staytime (30 days) is selected, and the 
optimum insertion and departure maneuvers (found by varying the orbit-insertion point for +90° to 
-90° true anomaly on the parking orbit) for that staytime are defined for each inclination. The 
insertion maneuver is then held fixed for each inclination, and the staytime is varied to determine 
the total AV requirement. Lrom these results, the available departure launch window for a 
given AV capability is determined. Reselection of the insertion maneuver would be expected to 
produce some gains, but based on the few cases investigated, changes that resulted in 
reduced AV at longer than nominal staytimes also resulted in increased AV for shorter than 
nominal staytimes. Such a tradeoff would be of interest during a study of a particular mission but 
adds little to the general discussion presented here. 

The minimum inclination that will allow coplanar arrival and departure is defined by the 
magnitude of the larger of the declinations at arrival and departure. This inclination is 28.7° and 
20.8° for the 1982 and 1986 missions, respectively. The maximum inclination that will allow 
coplanar arrival and departure is 180° minus the minimum inclination (i.e., a retrograde orbit). 
Therefore to obtain the parametric information contained in this report, the orbit inclination was 
varied from 30° to 1 50° for both missions in 10° intervals. Parking orbit eccentricities of 0.3, 0.5, 
and 0.7 were assumed with a periapsis altitude of 300 km. Both northern and southern insertion 
directions were considered. 

To decrease the volume of data but still outline the salient features of elliptical parking orbits, 
data are presented only for the better insertion direction of each year (i.e., northern insertion for 
1 982 and southern insertion for 1 986). In addition, neither data for an orbit eccentricity of 0.5 nor 
data for two-impulse transfers are presented since they resulted in no increased understanding of 
the problem and were generally bounded by the data shown. 
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Inclination, deg 


1982 inbound swingby — Figure 4(a) presents the total AV contours (arrival plus departure) 
for a one-impulse transfer from an elliptical orbit of 0.3 eccentricity as a function of inclination and 
staytime. If a coplanar periapsis arrival and departure were possible, a minimum AV of 6.6 km/sec 
could be achieved. As can be seen in figure 4(a), two small regions exist for which the arrival and 
departure conditions are sufficiently near the absolute optimum that essentially no penalty occurs. 
A 6-percent increase (0.4 km/sec) in the total AV results in a staytime of at least 20 days for all 
inclinations, and as much as 50 days for a polar orbit. Using the total AV representative of 
minimum launch windows for circular orbits (7.7 km/sec), a staytime of at least 40 days exists for 
all inclinations, compared with the 5-day staytime for the one-impulse transfer from circular orbits 
(ref. 2). 

Data for a three-impulse transfer from a 0.3 eccentricity orbit are shown in figure 4(b). Note 
that some increase in launch window duration for a given velocity capability occurs over that for a 
one-impulse transfer. However, the increase is not large due to the close positions of the arrival 
vector and the departure vector and their relative magnitudes. For this year, the arrival and 
departure pericenters are located so that both maneuvers can be performed with one impulse near 
the orbit pericenter because of the small plane-change requirement. A AV of 7.7 km/sec would 
produce a launch window of about 48 days for any inclination except those that are near 
equatorial. 



(a) 1982 inbound swingby; northern insertion; parking (b) 1982 inbound swingby; northern insertion; parking 
orbit eccentricity = 0.3; one impulse. orbit eccentricity = 0.3; three impulse. 

Figure 4.- Mars AV contours. 
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I 1 1 1 1 1 — I 1 I - I I 

o 20 40 60 80 100 

Staytime, days 

Figure 5.- Mars AV contours; 1982 inbound 
insertion; northern insertion; parking orbit 
eccentricity = 0.7; one impulse. 


Figure 5 presents the contour plot for a 
one-impulse transfer from a 0.7 eccentricity 
parking orbit. For this orbit, the minimum 
possible AV requirement is 5.45 km/sec, which is 
nearly achieved in two regions of figure 5. An 
increase of 5 percent in AV to 5.7 km/sec allows 
a window of 50 days for near polar and 
high-inclination posigrade orbits. The particular 
form of this AV contour in the retrograde region 
is of interest in that it indicates a maximum 
staytime of up to 50 days for a AV of 
5.7 km/sec, but also a minimum staytime of 20 to 
40 days for the same velocity capability. If this 
early departure restriction is of concern, then 
the AV requirements for retrograde orbits are 
higher than those for posigrade orbits for the same 
launch window. Use of a AV of 7.7 km/sec 
results in a launch window of 75 days for nearly 
all inclinations. 

The use of a three-impulse transfer from this 
orbit eccentricity does not signficantly change the 
launch window contour, for reasons given above; 
therefore it is not shown. 


1986 outbound, swingby - Data for a one-impulse transfer from a 0.3 eccentricity parking 
orbit is shown in figure 6(a). A minimum AV requirement of 6.0 km/sec exists for the idealized 
coplanar, pericenter maneuvers; a minimum of 6.1 km/sec was actually achieved. Since the arrival 
and the departure asymptotes differ over 45° in right ascension, a AV of 8.3 km/sec must be 
allowed before a generally available launch window exists for most inclinations at all staytimes. As a 
result, the one-impulse case requires a serious compromising of the departure maneuver. The 
difference in right ascension also results in the contour figure becoming nonsymmetric about the 
polar orbit. This nonsymmetry is such that retrograde orbits provide the maximum launch window 
for a given AV capability as can be seen in all the curves for the 1 986 mission. 


Figure 6(b) presents the contour maps for the three-impulse transfer. In this case, as compared 
to 1 982, the large plane-change requirement allows for effective use of the three-impulse technique 
and thus results in a significant increase in the available launch window. The use of a 
total AV representative of minimum launch windows for circular orbits (7.3 km/sec) results in a 
launch window of over 60 days for almost all inclinations. 


The contour map for a one-impulse transfer from a 0.7 eccentricity orbit is shown in 
figure 7(a). For this orbit, the minimum AV requirement would be 4.85 km/sec. Such alow value 
was not achieved anywhere on the one-impulse transfer plot. A AV of 7.3 km/sec opens the region 
below 50° and above 120° inclination for staytimes of 100 days. A AV of 8.0 km/sec allows use 
of the entire contour map. 

The results of using a three-impulse transfer from the 0.7 eccentricity orbit are shown in 
figure 7(b). Here a AV of 6.3 km/sec opens essentially the entire map to use, while a AV of only 
5.3 km/sec would allow any retrograde orbit. 
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Inclination, deg 




(a) 1986 outbound swingby; southern insertion; parking 
orbit eccentricity = 0.3; one impulse. 


(b) 1986 outbound swingby; southern insertion; parking 
orbit eccentricity = 0.3; three impulse. 


Figure 6.- Mars AV contours. 


Inclination, deg 




(a) 1986 outbound swingby; southern insertion; parking (b) 1986 inbound insertion; northern insertion; parking 

orbit eccentricity = 0.7; one impulse. orbit eccentricity = 0.7; three impulse. 


Figure 7.- Mars AV contours. 


Launch Window Comparison 


The previous sections have indicated in a parametric manner the effects of orbital parameters, 
transfer techniques, and interplanetary trajectories on the available launch window (i.e., allowable 
staytime variation) at Mars. In this section, the data are compared with those previously published 
for circular orbits (ref. 2) to provide insight into the potential tradeoffs among the pertinent system 
parameters (i.e., orbit eccentricity, energy requirements, number of impulses, trajectory, and 
staytime). For the comparison, a fixed total propulsive velocity requirement (arrival plus departure) 
of 8.7 km/sec was used for the 1982 opportunity. Figure 8 shows that for circular orbits, use of 
single impulse at this energy level allows staytimes of up to 35 days for all inclinations and up to 
70 days for specific inclinations. Use of three impulses makes the 70-day staytime available at all 
inclinations. With an orbit of intermediate eccentricity (e = 0.3), all inclinations have allowable 
staytimes up to 80 days with the optimal number of impulses depending on the particular 
inclination. This overlap is due to the method of selecting the arrival maneuvers and the small 
difference between the right ascensions and declinations of the arrival and departure vectors for the 
1982 mission. Use of a highly eccentric orbit (e = 0.7) results in essentially identical performance 
for the one- and three-impulse transfers. Staytimes of 90 days are readily available. 

I (a ) e = 0 (b) e = 0.3 (c ) e = 0.7 

Declination 
limitation 


and 3 


Declination 
limitation 

LuJ 1 l 1 I 1 l l l 

O 20 40 60 80 60 80 100 80 100 

Staytime, days 



Figure 8.- Effect of eccentricity on launch window 1982 inbound swingby; total AV = 8.7 km/sec northern 

insertion. 

Comparison of the circular and highly elliptical orbits reveals that for specific values of 
inclination they have nearly the same staytimes if a single impulse is employed; however, only the 
highly elliptical orbit allows the use of all inclinations at this energy level for the maximum 
staytime. Use of three impulses allows circular orbits to provide nearly the same staytime variation 
as the one-impulse elliptical. The major advantage of elliptical orbits is that the 


12 


total AV requirement can be reduced below the circular orbit level for the same staytime. 
Referring back to figure 5, it can be seen that for e = 0.7, a AV of only 6.7 km/sec can produce 
an allowable staytime of 70 days for all inclinations. 


Thus, an elliptical orbit (e = 0.7) with a one-impulse transfer can provide the same launch 
window as a circular orbit with a three-impulse departure while allowing a savings of 2 km/sec in 
the propulsive velocity requirement. Another advantage is the reduced complexity afforded by a 
single velocity maneuver at arrival and one at departure, although the effect of the elliptical orbit 
on operations at the planet may be undesirable. 


The 1986 mission is representative of missions with significant variations between the arrival 
and departure vectors. For the 1986 mission, a total AV of 7.3 km/sec for arrival plus departure 
was used. Figure 9 compares the allowable staytime for circular and elliptical orbits at this energy 
level. The shaded areas represent those regions of orbit inclination and staytime where the sum of 
the arrival and departure velocity maneuvers requires a velocity capability greater than 7.3 km/sec. 
The contours are unsymmetrical because of the velocity vector differences discussed, and near polar 
orbits are not available at this energy level for the single-impulse transfers. There are staytimes 
exceeding 60 days, but these are not shown since the contour lines become complex when the 
central angle of the trajectory from Mars to Earth exceeds 1 80° (e.g., fig. 6). 

(a) e = 0 (b) e = 0.3 ( c) e = 0.7 



Staytime, days 



0 20 40 60 


Figure 9.- Effect of eccentricity of launch window 1986; AV = 7.3. 


The energy level of 7.3 km/sec does not allow a very wide range of staytimes for circular 
orbits using the single impulse transfer. A AV of 8.3 km/sec would allow use of most posigrade 
orbits (ref. 2). Also, use of the three-impulse transfer technique allows all but a few nearly polar 
retrograde orbits. Use of medium eccentricity (e = 0.3) opens the low-inclination region for 
single-impulse transfers at this energy level, while three impulses again allow use of all inclinations 
except an area near 30° inclination. Increasing the orbit eccentricity to 0.7 opens the regions of low 
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posigrade and almost all retrograde orbits to allow 60-day staytimes. Three-impulse transfers would 
open the entire region of inclination — staytime shown. As with the 1982 mission, use of the highly 
elliptical orbit allows lower energies to be considered. Referring to figure 7 reveals that almost the 
entire region is available for a AV of 6.3 km/sec and that all retrograde orbits are available for only 
5.3 km/sec where a three-impulse transfer is used. This is a significant reduction from the circular 
orbit requirements and even from the single-impulse elliptical orbit. 

The use of a three-impulse transfer for any orbit eccentricity provides flexibility in the 
allowable inclination for specific staytimes. This flexibility in choice of inclination has great 
implications for mission analysis, for it means that the orbit selection and the interplanetary 
trajectory selection can be decoupled. The orbit therefore can be selected on the basis of scientific 
goals, experiment requirements, and/or operational requirements without being compromised by 
interplanetary trajectory considerations. 


CONCLUDING REMARKS 


Increasing eccentricity of elliptical orbits reduces the total AV requirement for a given 
launch window (staytime). In addition, the use of three impulses to perform the transfer generally 
requires lower AV than does the use of a single impulse. However, orbits can exist (e.g., 1982, 
e = 0.7) for which no gain is shown by the three-impulse maneuver. This situation usually 
corresponds to a low plane-change requirement. Two specific conclusions can be drawn. First, as a 
result of the effectiveness of the three-impulse transfer in reducing the AV requirement to a nearly 
constant value regardless of the orbit geometry, the orbit selection can be optimized from other 
considerations (i.e., science, communication, power, etc.) independent of the interplanetary 
trajectory considerations. Second, for highly elliptical orbits (e = 0.7) a total AV (arrival plus 
departure) at Mars that is 1 km/sec below the minimum coplanar requirement for circular parking 
orbits provides staytimes of 60 days and usually allows any orbit inclination to be used. This AV is 
effectively 2 km/sec below that required for circular orbits with the same allowable staytime. In 
addition, whereas the circular orbits require three impulses at departure, for some missions a single 
impulse can be used to eject from the elliptical orbit. 


National Aeronautics and Space Administration 

Moffett Field, California 94035, May 14, 1970 
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APPENDIX A 


SINGLE IMPULSE ANALYSIS 


Figure 10 indicates the geometry and notation used in the analysis of single impulse escape 
maneuvers. At the time of departure, the elliptical orbit plane is specified by the longitude of the 
ascending node £2 (measured from planet Vernal Equinox) and the inclination i. Periapsis of the 
orbit in this plane is specified by the argument of periapsis co measured in the orbit plane from the 
ascending node, and is designated by the unit vector r p . The departure asymptote is specified by 
two angles: the right ascension p and declination 8 . The pertinent unit vectors n (the normal to 
the orbit plane) and V„ are given by 

n = sin ft sin i i - cos ft sin i j + cos i k 
V = cos p cos 6 i + sin p cos 6 j + sin 6 k 

CO 

and the unit vector r p is obtained from 

r n x x p = sin u n 
r n ' r p = cos m 

where r n is the unit vector at the ascending node and is given by 

x n = cos ft i + sin ft j 


z 
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Given a true anomaly of departure d (see 
fig. 1 1 ), the radius of the escape hyperbola is 
equal to the radius of the ellipse at that point. 



r 


a 0 Cl - e Q 2 ) a^l - e h 2 ) 

1 + cos 0 ~ 1 + e, cos 0. 
0 h h 


(Al) 


where 


a h = - 



The angle $ between the radius vector, r, and 
the V«, vector is given by 

cos <|> = r • V OT (A2) 

where f is determined from 


Figure 1 1 .- Ejection point geometry. 


and 


r x r = sin 0 n 
P 


r • r = cos 0 
P 


Equation (A2) is not explicit, since there are two transfer possibilities for a given point on an ellipse 
to a given hyperbolic asymptote. One involves a rotation from r to V„ in a positive sense, and the 
other involves such a rotation in a negative sense. Thus, the angle <j> determined from 
equation (A2) has two values, one plus and one minus. The orientation of the two possible 
departure hyperbolas is shown in figures 12(a) and (b). 


The true anomaly on the hyperbola for departure maneuver 6\y is then given by 

0 = it - <j> - e 

h 


(A3) 


where 


e 


cos -1 1/e. 

h 


Noting from equation (A3) that 


cos 0^ = - cos (cj> + e) 


cos <j> cos e + sin <j> sin e 


cos <j> 


J\ - e^ 2 sin <j> 
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(b) Negative rotation. 

Figure 12.- Orientation of departure hyperbola. 


and substituting into equation (A1 ) yields 

V 1 ~ e h 2) 

1 - cos <t> + J- - e, 2 sin ((> 

This equation can be solved for the escape hyperbola eccentricity ejj by substituting 
? = >/ e h 2 ~ 1 ar *d solving the resultant quadratic equation for £. 

-r sin (j) ± J r 2 sin 2 (J> - 4a^(l - cos <j)) 
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The positive radical is not a valid solution since the semimajor axis for a hyperbola is negative and f 
is by definition positive. Therefore, 


e h = 


1 - 


r sin <j> 
2a h 


^/r 2 sin 2 <j) - 4a^(l - cos 40 


1/2 


2a h 


The velocity increment vector direction and magnitude to change from the parking elliptical 
orbit to the escape hyperbola can now be determined through vector analysis as follows. The 
direction cosines of the unit velocity vector on the ellipse at the point of departure 6 are 
determined by 


r x V Q = sin - y 0 )n 
and 

* • V Q = cos ( j - y 0 ) 


where the flight path angle y Q is given by 


e sin 0 

tan y - j f e C os 6 
o 

Similarly, the direction cosines of the desired unit velocity vector on the escape hyperbola can be 
determined after the true anomaly on the departure hyperbola is determined from equation (A3) 
and the unit vector normal to the departure transfer plane is determined from 

r x V 

/v 00 

n, = : 

h sin <j> 

Then, the unit velocity vector on the hyperbola is given by 

r X V h = s l n (j - T h ) n h 

and 

i ■ \ ‘ cos(l - Y h ) 
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where the flight path angle 7^ is given by 

tan Yi_ 
h 


e, sin 0, 
h h 


1 + e. cos 0, 
h h 


The required velocity increment is then given by 


where 


and 


AV 


= V, V, - 
h h 


V V 
o o 






Since the transfer in a negative rotation sense from a given point on an ellipse to a given hyperholic 
asymptote involves passing through periapsis on the departure hyperbola, it is necessary to calculate 
the periapsis altitude to determine if a positive altitude exists for this transfer. The periapsis altitude 
is given by 

H p = V 1 ‘ e h) " R 

where R is the radius of the planet. 

Finally the optimum departure position on the ellipse is determined by numerical search by 
stepping true anomaly around the parking ellipse. 
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APPENDIX B 


TWO-IMPULSE ANALYSIS 


Figure 10 indicates the geometry and notation used in the analysis of the restricted 
two-impulse escape maneuver. The description of the geometry at departure is the same as the 
single-impulse analysis contained in appendix A. The restricted two-impulse maneuver studied here 
involves an impulse at apoapsis to rotate the parking ellipse coplanar with the departure asymptote 
followed by a tangential impulse to transfer in that plane from the parking ellipse to the hyperbolic 
asymptote. 

The departure point on the ellipse as shown in figure 13 is designated by the true 
anomaly 0 - L and also by the true anomaly on the escape hyperbola 0^. The key independent 
parameter is the angle between the periapsis unit vector ?p and the departure asymptote unit 
vector V OT and is designated as P. It is convenient to invert the escape problem and the roles of 
dependent and independent parameters by specifying 0j and then determine the resulting angle P. 



Figure 13.- Escape geometry. 
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Applying the conditions of tangency between the escape hyperbola and the intermediate 
ellipse at the point where AV is applied and then solving for the unknown eccentricity of the 
hyperbola yields: 


e 


h 


1 



1/2 


where 


r = 


r (1 + e-) 

P 1 

1 + cos 0^ 


\ = - 


P 



and 


sin 0^ 

tan y = -z 

1 + cos 


The true anomaly on the escape hyperbola is then given by 

13 

'h r/a h 


1 - e 2 

h tan y 

sm 0, = — — - 

h e. 


COS 0, = — 

h e. 


1 - V 

r/a h 


- 1 


Referring to figure 13, the asymptote half angle is determined by 

-1 1 
e = cos — 

e h 

and the angle P is finally determined by 

P = 7T-£+0. - 0 
i h 

The second velocity increment AV 2 is given directly by 


AV, 


2 „ 2 
U — + V 

X co 


1/2 


1 - e. 


1/2 
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With reference to figure 10, the unit vector normal to the original orbit plane n is given by 

A A A A 

n = sin 0 sin i i - cos 0 sin i j + cos i k 

The unit vector, rp, is then determined by 

r x r = sin am 
n p 


and 


r • r = cos a) 
n p 

where f n is the unit vector at the ascending node and is given by 

r n = cos 0 i + sin 0 j 

The unit vector in the direction of the departure asymptote is given by 

V = cos 6 cos p i + cos 6 sin p j + sin 5 k 

OO ^ J 


. and the angle P is determined from 


cos P = r • V 

p 00 

The unit vector normal to the departure orbit plane (see fig. 1 3) is then given by 



r x V 

_P 


00 


sin P 


and thus the plane change angle required at apoapsis to rotate the original orbit plane coplanar 
with is given by 

A A 

cos A = n^ • n 

The required velocity increment is 


AV 


1 




(1 - e i) 


1/2 


sin 


The total velocity requirement for the escape maneuver is thus given by 

AV = AV X + AV 2 
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